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Charmless 3-body decays of B mesons are studied in a simple model based on the factorization 
approach. We have identified two distinct sources of nonresonant contributions: a small contribution 
from the tree transition and a large source of the nonresonant signal in the matrix elements of scalar 
densities, e.g., (KK\ss\0), induced from the penguin transition. This explains the dominance of the 
nonresonant background in B —> KKK decays, the sizable nonresonant fraction in K~-k + tt~ and 
K it + tv~ modes and the smallness of nonresonant rates in B — > tttttt decays. The seemingly huge 
discrepancy between BaBar and Belle for the nonresonant contribution in the decay B~ — > K~tt tt~ 
is now relieved. We have computed the resonant and nonresonant contributions to charmless 3- 
body decays and determined the rates for the quasi-two-body decays B — > VP and B — > SP. 
Time-dependent CP asymmetries sin 2/3 e ff and Acp in K + K~ Ks, KsKsKs, Kstt + tt~ and i\s7r 7r° 
modes are estimated. 



I. INTRODUCTION 

Recently many three-body B decay modes, for ex- 
ample, K+tt+tt'-, K°7r+7r-, K+tt~tt°, K+K+R-, 
K°K + K~ , KsKsK + and tt + tt + tt~ , have been ob- 
served at B factories with branching ratios of order 
10~ 5 . The Dalitz plot analysis of 3-body B decays 
provides a nice methodology for extracting informa- 
tion on the unitarity triangle in the standard model. 
The three-body meson decays are generally dominated 
by intermediate vector and scalar resonances, namely, 
they proceed via quasi-two-body decays containing a 
resonance state and a pseudoscalar meson. Indeed, 
most of the quasi-two B decays are extracted from 
the analysis of three-body B decays using the Dalitz 
plot technique. Three-body hadronic B decays in- 
volving a vector meson or charmed meson also have 
been observed at B factories. In this talk I'll focus on 
charmless 3-body B decays. 

Experimentally, there are two striking features for 
3-body hadronic B decays: 

(i) large noresonant fractions in peguin-dominated 
modes 

It is known that the nonresonant signal in charm de- 
cays is small, less than 10% I n the past few years, 



TABLE I: Fractions (%) of nonresonant contributions to 
various charmless three-body decays of B mesons It 
will be explained below about the BaBar measurement of 
the nonresonant component in the B~ — > K~tt + tt~ decay. 



Decay 


BaBar 


Belle 


B~ -> K + K~K~ 


141 ± 18 


74.8 ± 3.6 


B° -» K+K-K° 


112 ± 15 




B° -» TTVtt- 




41.9 ± 5.5 


B~ -> K-n+iT- 


17.lt™ 


34.0 ±2.9 


B° -» K+7V-7T 




15.6 ± 7.7 < 25.7 


B~ —* n + n~n~ 


13.6 ±6.1 





some of the charmless B to 3-body decay modes have 
been measured at B factories and studied using the 
Dalitz plot analysis. We see from Table U that the 
nonresonant fraction is about ~ 90% in B — > KKK 
decays, ~ 17 — 40% in B — > Kitty decays (smaller in 
the Ktttt decay), and ~ 14% in the B — > iririr de- 
cay. Hence, the nonresonant 3-body decays play an 
essential role in penguin-dominated B decays. While 
this is a surprise in view of the rather small nonreso- 
nant contributions in 3-body charm decays, it is not 
entirely unexpected because the energy release scale 
in weak B decays is of order 5 GeV, whereas the ma- 
jor resonances lie in the energy region of 0.77 to 1.6 
GeV. Consequently, it is likely that 3-body B decays 
may receive sizable nonresonant contributions. It is 
important to understand and identify the underlying 
mechanism for nonresonant decays. 

Nonresonant amplitudes in charm decays are usu- 
ally assumed to be uniform in phase space. However, 
this is no longer true in B decays due to the large 
energy release in weak B decays. While both BaBar 
and Belle have adopted the paramctrization 

i NR = (ci 2 e i * 12 e- Qfl " ±ci 3 e* 3 e- QS i3 

+ c 23 e l ^e~ as ^)(l + b NR e^ +s ^) (1) 

to describe the non-resonant B — » KKK amplitudes, 
they differ in the analysis of the nonresonant com- 
ponent in B — > Ktttt decays. Belle still employed 
the above exponential parametrization to analyze the 
nonresonant contribution, but BaBar used the LASS 
parametrization to describe the Kit 5-wave and the 
nonresonant component by a single amplitude sug- 
gested by the LASS collaboration to describe the 
scalar amplitude in elastic Ktt scattering. Since the 
BaBar and Belle definitions of the .Kq(1430) and non- 
resonant differ, the branching fractions and phases are 
not directly comparable. We will come this point more 
in Sec. III. 

Experimentally, it is hard to measure the direct 3- 
body decays as the interference between nonresonant 
and quasi-two-body amplitudes makes it difficult to 
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disentangle these two distinct contributions and ex- 
tract the nonresonant one. 

(ii) New broad scalar resonances /x(1550) and 
/x(1300) 

A broad scalar resonance /x(1500) (or Vo(1550) 
denoted by BaBar) has been seen in B — » K + K + K~ , 
K + K~Ks and K + K~tt + decays at energies around 
1.5 GeV. However, it cannot be identified with the 
well known scaler meson Jo (1500). This is because 
/o(1500) decays into 7r + 7r _ about five times more fre- 
quently than to K + K~ . Identification of /x(1500) 
with /o(1500) will imply that the K + K~ peak at 
1.5 GeV will be accompanied by a peak in 7r + 7r~, 
which is not seen experimentally. Hence, the nature 
of /x(1500) is not clear. 

Moreover, there exists a production puzzle for 
/x(1500). Both BaBar and Belle have seen a 
large fraction from /x(1500) in the decay B~ — > 
K+K-R-; (121 ±19 ±6)% by BaBar and (63.4 ± 
6.9)% by Belle [J], whereas it is only about 4% seen 
by BaBar in B° -> K+R-Rs [|. The puzzle is 
that why /x(1500) behaves so dramatically different 
in charged and neutral B decays to 3 kaons. It is not 
clear whether the large production of /x(1500) is a 
genuine effect or just a statistical fluctuation. Any- 
way, this issue should be clarified soon. Notice that 
Belle actually found two solutions for the fraction of 
ix(1500)if J in B~ -» K + K+K~ Q: (63.4 ± 6.9)% 
and (8.21 ± 1.94)%. The first solution is preferred by 
Belle. It is probably worth of re-examining the small 
solution. 



II. THREE-BODY B DECAYS 

In analog to two-body decays of heavy mesons 
which can be analyzed using the model-independent 
quark diagrammatic approach, three-body decays of 
the heavy mesons can be expressed in terms of some 
quark-graph amplitudes (see Fig. [I]): T\ and 72, 
the color-allowed external JF-emission tree diagrams; 
Ci and C2, the color-suppressed internal VF-emission 
diagrams; £, the W-exchange diagram; A, the W- 
annihilation diagram; Vi and V2, the penguin dia- 
grams, and V a , the penguin- induced annihilation dia- 
gram. The quark-graph amplitudes of various 3-body 
B decays B — > Tih + h~ and B — * Kh + h~ are summa- 
rized in Table I of 8] . As mentioned in Q , the use of 
the quark-diagram amplitudes for three-body decays 
are in general momentum dependent. This means that 
unless its momentum dependence is known, the quark- 
diagram amplitudes of direct 3-body decays cannot be 
extracted from experiment without making further as- 
sumptions. Moreover, the momentum dependence of 
each quark-diagram amplitude varies from channel to 
channel. 

We take the decay B — * K + K~K as an illustra- 
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FIG. 1: Quark diagrams for the three-body decays of 
heavy mesons, where Q denotes a heavy quark. 

tion. Under the factorization approach, its decay am- 
plitude consists of three distinct factorizable terms: 
(i) the transition process induced by b — > s pen- 
guins, (B° — > K + T?°) x (0 — > K~), (ii) the current- 
induced process through the tree b — > u transition, 
(75° -> if) x (0 -> K + K~), and (iii) the annihilation 

process (B° -> 0) x (0 -> K + K~~K°), where (A -> B) 
denotes a A —> B transition matrix element. 



A. Nonresonant background 

For the transition process, the general expression of 
the nonresonant contribution has the form 

{K-(p 3 )\(su) v ^ A \0){K (pAK+(p 2 )\(ub) v _ A \B ) NR 

= \2m\r + (m| - s i2 - m 2 Auj + 

+ (s 2 3 ~ sia)w-] , (2) 

where (qiq 2 )v-A = 9i7c(l ~ 75)92- In principle, one 
can apply heavy meson chiral perturbation theory 
(HMChPT) to evaluate the form factors r, uj + and 
uj- (for previous studies, see @). However, this will 
lead to too large decay rates in disagreement with 
experiment [Tfjj . A direct calculation indicates that 

the branching ratio of B — > K + K~K arising from 
the transition process alone is already at the level of 
77 x 10 -6 which exceeds the measured total branching 
ratio of 25 x 10~ 6 . The issue has to do with the 
applicability of HMChPT. In order to apply this ap- 
proach, two of the final-state pseudoscalars (K + and 

K° in this example) have to be soft. The momentum 
of the soft pseudoscalar should be smaller than the 
chiral symmetry breaking scale A x of order 0.83 — 1.0 
GeV. For 3-body charmless B decays, the available 
phase space where chiral perturbation theory is appli- 
cable is only a small fraction of the whole Dalitz plot. 
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Therefore, it is not justified to apply chiral and heavy 
quark symmetries to a certain kinematic region and 
then generalize it to the region beyond its validity. If 
the soft meson result is assumed to be the same in 
the whole Dalitz plot, the decay rate will be greatly 
overestimated. 

Recently we have proposed to parametrize the b — > 
u trasnition-induced nonresonant amplitude given by 
Eq. © as [HI 



A 



NR 



.4 



HMChPT -a. 
NR e 



■Oi+p 2 ) e i0i 



(3) 



so that the HMChPT results are recovered in the 
chiral limit p\, pi — > 0. That is, the nonreso- 
nant amplitude in the soft meson region is described 
by HMChPT, but its energy dependence beyond 
the chiral limit is governed by the exponential term 
e - a NR PB-{pi+P2) _ -jijjg un k n own parameter a NR can be 
determined from the data of the tree-dominated decay 
B~ — > 7r + 7r~7T~. 

In addition to the b — > u tree transition, we need to 
consider the nonresonant contributions to the b — > s 
penguin amplitude 



A-2 



{K° | (sb) V -A\B°) (K + K- 1 (uu) V -a\0) , 
(K°\sb\B°)(K + K-\ss\Q). 



(4) 



away from the resonant region. In the presence of the 
resonances, we thus need to introduce a nonresonant 
<t nb term which can be constrained by the measured 

pf — > KsKsKs rate and the K + K~ mass spectrum 



B. Resonant contributions 

Vector meson and scalar resonances contribute 
to the two-body matrix elements (-Pif^l^O) and 
(-P1-P2ISIO), respectively. They can also contribute 
to the three-body matrix element {P\P2\V ll — A^B). 
Resonant effects are described in terms of the usual 
Brcit-Wigner formalism. More precisely, 

(K+K-\q^ q \Q) R = J2(K + K-\Vi) 

i 

1 



TOy. — s — iniVi 



(K+K-\ss\0) R = ^2{K + K-\Si) 

i 

1 



m 2 s . — s — imsj^s., 



■{Si\ss\0}, 



(8) 



The 2-kaon creation matrix elements can be expressed 
in terms of time-like kaon current form factors as 



(K+{p K+ )K-{p K -)\q^q\0) = (p K+ -p K -) 
(K O (p K o)K a (p R o)\qjf,q\0) = {pk«-Pk<> 



F K+K- 



where Vi = (f>, p, u>, ■■■ and Si = 
/o(980),/o(1370),/ (1500),---. In this manner 
we are able to figure out the relevant resonances 
which contribute to the 3-body decays of interest and 
compute the rates of B — > VP and B — > SP. 



The weak vector form factors F^ +K and F^° K ° can 
be related to the kaon electromagnetic (e.m.) form 
factors F^ K and F^°^° for the charged and neu- 
tral kaons, respectively. Phenomenologically, the e.m. 
form factors receive resonant and nonresonant contri- 
butions 



f k+k- _ F 

± em ± p 



F u + Fa, 



F 



K U K U 



-F„ + F ul 



F,, 



NR, 

- f 'nr- 



(6) 



The resonant and nonresonant terms in Eq. © can 
be determined from a fit to the kaon e.m. data. 
The non-resonant contribution to the matrix element 
(K + K~ \ss\0) is given by 



(K+(p 2 )K-(p 3 )\ss\0) 



NR 



^F'nr) 



(7) 



The nonresonant cr NR term is introduced for the fol- 
lowing reason. Although the nonresonant contribu- 
tions to ff- K and Fg K are related through the equa- 
tion of motion, the resonant ones are different and 
not related a priori. As stressed in [12j], to apply 
the equation of motion, the form factors should be 



III. PENGUIN-DOMINATED B 
B -> Ktvtt DECAYS 



KKK AND 



As mentioned in the previous section, we employ 

the decays B — > K + K~K and KsKsKs to fix the 
nonresonant parameter er NR to be 



= e 



,171-74 



(3.36±S;£) GeV. 



(9) 



It turns out that the nonresonant contribution arises 
dominantly from the transition process (88%) via 
the scalar-density-induced vacuum to K K transition, 
namely, (K + K~\ss\0), and slightly from the current- 
induced process (3%). Physically, this is because the 
decay B — > KKK is dominated by the b — > s penguin 
transition. The nonresonant background in B —> KK 
transition does not suffice to account for the experi- 
mental observation that the penguin-dominated decay 
B — > KKK is dominated by the nonresonant contri- 
butions. This implies that the two-body matrix el- 
ement e.g. (KK\ss\0) induced by the scalar density 
should have a large nonresonant component. 

We have considered other B —> KKK decays such 
as B~ -> K+K-K- and B~ -> K-K S K S and 
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TABLE II: Branching ratios (in units of 10 6 ) of resonant and nonresonant (NR) contributions to B — ► K tv + tv . 
Theoretical errors correspond to the uncertainties in (i) a NR , (ii) m s , F SK and a NR , and (iii) 7 = (59 ± 7)°. For 

the BaBar results, the branching fraction of K (1430)7r~ comes only from the Breit-Wigner component of the LASS 
parametrization, while the nonresonant contribution includes both the nonresonant part of the LASS shape and the 
phase-space nonresonant piece. 



Decay mode 


BaBar [13] 


Belle [14] 


Theory [11] 


K*\- 


7.2±0.4±0.7±£'j 


6.45 ± 0.43 ± 0.48^;^ 


n +o.o + o.s + o.o 

°- u -0. 0-0. 7-0.0 


7?o°(1430)7r- 


32.0 ± 1.2 ± 2.7±Ia ± 5-2 


32.0 ± 1.0 ± 2.4jl^| 


m c+o. 0+3. 2+0.0 

—0.0 — 2.7—0.1 


p°K- 


3.56 ± 0.45 ± 0.43to;?s 


3.89 ±0.47 ±0.291^9 


1 O+0. 0+1. 9 + 0.1 

1 - J -0. 0-0. 7-0.1 


fo(9SO)K- 


10.3 ± 0.5 ±i.3tl:l 


8.78 ± 0.82 ± 0.65^1:11 


7 7+0.0+0.4+0.1 
'■'-0.0-0.8-0.1 


NR 


9.3±1.0±1.2±o;4±1.2 


16.9 ± 1.3 ± 1.3±o.g 


1 7+0.5+11.0+0.2 
10 - ' -0.6- 6.3-0.2 


Total 


54.4 ± 1.1 ±4.6 


48.8 ±1.1 ±3.6 


a r n +0. 3+16. 4+0.1 
<±J.U_ 4 _ 10 s _o j 



TABLE III: Branching ratios (in units of 10 6 ) of resonant and nonresonant (NR) contributions to B —* K -k + tv 



Decay mode Belle [16] Theory [11] 

K*-* + 5.6 ± 0.7 ± o.5±S:| 2.iini ; 5 5 t : 3 3 

^-(1430)rr+ 30.8 ± 2.4 ± 2.4±°; 8 lO-HH^oi 



o°K u 6.1±1.0±0.5±i° 2 0+0.0+1.9+0.1 



/o(980)F° 7.6±1.7±0.7lg; 5 7 7 .7 + . - + - A 7 + - 

NR 19.9 ±2.5 ±1.61^ 15-6tg;£i;|tg;g 

Total 47.5 ± 2.4 ± 3.7 ^.0t°'^t lol- 0.0 



found that they are also dominated by the nonres- 
onant contributions. Our predicted branchin g; ra tio 
B(B- -» K+K-K-) NR = (25.3+^) x 10" 6 [ll| is 
in good agreement with the Belle measurement of 
(24.0^1°) x 10_6 > but a factor of 2 smaller than the 
BaBar result of (50 ± 6 ± 4) x 1CT 6 0. 

The resonant and nonresonant contributions to the 
decay B~ — ► K~n + ir~ arc shown in Table HU We 
see that the calculated K *ir and pK rates are smaller 
than the data by a factor of 2 ~ 3. This seems to be 
a generic feature of the factorization approach such 
as QCD factorization where the predicted penguin- 
dominated VP rates are too small compared to ex- 
periment. We shall return back to this point later. 

At first sight, it appears that the nonresonant 
branching ratio (2.4 ± 0.5 ± l-3±g;l) x 10~ 6 in B~ 
i^~7r + 7r _ obtained by BaBar [13j is much smaller than 
the one (16.9± 1.3± 1.3^g) x 10" 6 measured by Belle 
[lij . However as mentioned in the Introduction, since 
the BaBar and Belle definitions of the A'q(1430) and 
nonresonant differ, it does not make sense to compare 
the branching fractions and phases directly. While 
Belle [lij employed the exponential parametrization 
Eq. Cp t o describe the nonresonant contribution, 
BaBar [13| used the LASS parametrization to describe 
the K 7T 5-wave and the nonresonant component by a 
single amplitude suggested by the LASS collaboration 



m r ^ 



M = 



m KT[ 



(m, 



l Kt. 



im T 



(10) 



q cot Si 



><1 



where cot 8b = ~^ + \ r Q- Since the LASS parametriza- 
tion is valid (experimentally confirmed) up to the Kit 
invariant mass of order 1.8 GeV, BaBar introduced 
a phase-space nonresonant component to describe an 
excess of signal events at higher Kir invariant mass. 
Hence, the BaBar definition for the ^(1430) includes 
an effective range term to account for the low Kn S- 
wave while for the Belle parameterization, this com- 
ponent is absorbed into the nonresonant piece. To 
stress once again, the result B(B~ — > K~tt + 7t~) n:r = 
(2.4 ± 0.5 ± 1.3±g:|) x 10~ 6 cited by BaBar is solely 
due to the phase-space nonresonant piece. 

From the above discussion, it is clear that part 
of the LASS shape is really nonresonant which has 
a substantial mixing with Kq(1A30). In principle, 
this should be added to the phase-space nonreso- 
nant piece to get the total nonresonant contribu- 
tion. Once this is done, it is possible that BaBar 
and Belle might agree with each other. Indeed, very 
recently BaBar have carried out this task [l3j]. By 
combining coherently the nonresonant part of the 
LASS parametrization and the phase-space nonreso- 
nant, BaBar found the total nonresonant branching 
fraction to be (9.3 ± 1.0 ± 1.2±®;| ± 1.2) x 10" 6 with 
the fit fraction being (17.1 ± 1.7 ± 1.6± x £f)% 
We see from Table |H] that the BaBar result is now 
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6 — o 
TABLE IV: Branching ratios (in units of 10 ) of resonant and nonresonant (NR) contributions to B 

the BaBar measurement, the resonance #o(1430) is replaced by the S-wave Ktt state, namely, (Ktt)q. 


-> if-7r + 7r . For 


Decay mode 


BaBar [18] 


Belle [19] 


Theory [11] 


K TT 

JsT*-(1430)7r + 
Z;°(1430)7r° 
p+K~ 
NR 


4.2+^5 ±0.3 
2.4 ±0.5 ±0.3 
9AtHtli ±1.8 

o 7+I.I + I.8 i o o 
°- ' -0.9-1.3 - 1 - 

s.ot°i, ± °- 6 


A q+1.5+0.5+0.8 
^•^-1. 5-0. 3-0. 3 

< 2.3 

5.1 ± l-5±g;? 

a 1 +1.6+0.5 
"•^-l. 5-0.6 
1 e: 1+3.4+1.4+2.0 
xo - 1 -3.3-l. 5-2.1 
5 7 +2.7+0.5 Q 4 
J - ' -2.5-0.4 *- 


i n+0. 0+0. 3+0.1 
1 - u -0. 0-0. 3-0.1 
i n +0.0+0.3+0.2 
i,U -0. 0-0. 2-0.1 
r n +0. 0+1. 5+0.1 
"•"-0.0-1.3-0.1 
A 9+O.O+I.4+O.O 
^■^-0. 0-1. 2-0.0 
9 r+0.0+3.6+0.2 
z '°-0.0-1.4-0.2 
n a+0. 3+6. 6+0.0 
y - u -0. 2-3. 5-0.0 


Total 


35-7±?-J ± 2.2 


36.6l4 i ±3-0 


90 q+0. 2+16. 1+0. 2 
ZO.a_ 2 _ 9.4-0.2 



TABLE V: Same as Table |n] except for B~ — > Tr + 7r Tf . The nonresonant background is used as an input to fix the 
parameter q nr defined in Eq. Q. 



Decay mode 


BaBar [20] 


Theory [11] 


P °1T- 


8.8 ± 1.0 ±0.6^7 


7 7+0.0+1.7+0.3 
'•'-0.0-1.6-0.2 


/o(980)rr- 


1.2 ±0.6 ±0.1 ±0.4 < 3.0 


r, oQ+0. 00+0. 01+0. 03 
u -" 3 -0. 00-0. 01-0. 02 


NR 


2.3 ±0.9 ±0.3 ±0.4 < 4.6 


input 


Total 


16.2 ± 1.2 ±0.9 


i 9 n +l.l+2.0+0.4 
1Z,.U_ 1 .2-1.8-0.3 



consistent with Belle within errors, though the agree- 
ment is not perfect as BaBar and Belle have different 
models for the nonresonant Kir mass. Likewise, the 
BaBar branching fraction (24.5 ± 0.9 ±2.1+£$) x 10" 6 
for B~ -> KS(U30)tt- cited in Table II of [3 in- 
cludes an effective range nonresonant component. In 
order to compare with the Belle result determined 
from the Breit-Wigner parametrization, it would be 
more appropriate to consider the Breit-Wigner com- 
ponent only of the LASS parametrizaion. The result 
is B(B~ -> #£(1430)71--) = (32.0 ± 1.2 ± 2.7±l\\ ± 
5.2) x 10 -6 which is now in good agreement with the 
Belle measurement (see Table |LT| . 

From Table [TT] we see that our predicted nonreso- 
nant rates are consistent with the Belle and BaBar 
measurements within errors. The reason for the 
large nonresonant rates in the K~n + ir~ mode is 
that under SU(3) flavor symmetry, we have the re- 
lation (Kn\sq\0) NR = (KK\ss\0) NR . Hence, the 
nonresonant rates in the A"~7r±7r~ (Table HI)) and 

K ir + Tr~ (Table HLT)) modes should be similar to that 

in K+K-K° or K + K~K~. Since the KKK chan- 
nel receives resonant contributions only from <fi and /o 
mesons, while K* , Kq, p, f resonances contribute to 
Kirn modes, this explains why the nonresonant frac- 
tion is of order 90% in the former and becomes of order 
40% in the latter. It is interesting to notice that, based 
on a simple fragmentation model and SU(3) symme- 
try, Gronau and Rosner [17j also found a large non- 
resonant background in K~'k + tt~ and K it + tt~ . 

Recently, BaBar has reported a new Dalitz-plot 
analysis of the decay B — ► K~tt + -k [l8j] (see Table 



IIV|) . Just as the K tt + tt mode, the reported non- 
resonant branching fraction B{B — > K~ir + ir°) nl = 
(4.4 ± 0.9 ± 0.5) x 10~ 6 by BaBar is only the phase- 
space part of nonresonant contributions. To get the 
total nonresonant rate, it is necessary to add the 
nonresonant component of the LASS parametriza- 
tion to the phase-space piece. When this is done, it 
will be interesting to compare the measured nonreso- 
nant branching fraction with our prediction B{B — > 

#-7r + 7r°) m . = (Q.elgjt!;!!) x 10 ~ 6 - lt should be 

stressed that the measured partial rates for B° — > 
(Kir)*-TT+ and (Kn)* n by BaBar (see Table|jyj in- 
clude an effective range Kit nonresonant component. 
Hence, it is not pertinent to compared them directly 
with the respective Belle measurements. 



IV. TREE-DOMINATED B ivnn, KKn 
MODES 

The B — > 7r7r7r mode receives nonresonant contri- 
butions mostly from the b — ► u transition as the non- 
resonant contribution in the penguin matrix element 
(7T + 7r _ |dd |0) is suppressed by the smallness of penguin 
Wilson coefficients a§ and a$. Hence, the measure- 
ment of the nonresonant contribution in this decay 
can be used to constrain the nonresonant parameter 
a NR in Eq. ©. 

Note that while B~ — > 7r + 7r~7r~ is dominated by 
the p° pole (Table IV)) . the decay B — > 7r + 7r~7r° re- 
ceives p^ and p° contributions. As a consequence, the 
7T±7r~7r mode has a rate larger than 7r + 7r~7r~ even 
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though the former involves a ir° in the final state. We 
predict that B{B° -> tt+tt-tt ) « 26 x 10" 6 

Among the 3-body decays we have studied, the de- 
cay B~ — > A + A~7r~ dominated by b — > u tree tran- 
sition and b — > d penguin transition has the small- 
est branching ratio of order 4 x 10~ 6 . BaBar [2l[ 
has recently reported the observation of the decay 
B + — > A + A~7r + with the branching ratio (5.0 ± 
0.5 ± 0.5) x 10 -6 . Our prediction for this mode, 
(4.0to'g±o sio 3) x 10 ~ 6 > is in accordance with exper- 
iment. 

V. QUASI-TWO-BODY B DECAYS 

It is known that in the narrow width approximation, 
the 3-body decay rate obeys the factorization relation 

T{B -> RP -> P1P2P) = T{B -> RP)B(R -> P1P2), 

(11) 

with R being a vector meson or a scalar resonance. 
Using the experimental information on B(R — > 

a 



£(A*° -> A+tT) 

B(A*°(1430) -> X+7T-) 
£(0 -> A+A") 



= £(A* + A°7T + ) 

= 2£(A* + -> A + 7T°) = |, 

= 2S(A*+(1430) -> A+tt ) 

= ^(0.93 ±0.10), 

= 0.492 ± 0.006, (12) 



we have extracted the branching ratios of B — > VP 
and P — > SP. The results are summarized in Table 
ED The predicted pir, / (980)A and / (980)tt rates 
are in agreement with the data, while the calculated 
<j>K, K*w, pA and Aq(1430)7t are in general too small 
compared to experiment. The fact that this work and 
QCDF lead to too small rates for </>A, K*ir, pK and 
Aq (1430)7r may imply the importance of power correc- 
tions due to the non-vanishing pa and p^ parameters 
arising from weak annihilation and hard spectator in- 
teractions, respectively, which are used to parametrize 
the endpoint divergences, or due to possible final-state 
rescattering effects from charm intermediate states 
22]. However, this is beyond the scope of the present 
work. 



VI. TIME-DEPENDENT CP ASYMMETRIES 

The penguin-induced three-body decays P° — > 
K + K~Ks and A5A5A5 deserve special attention as 
the current measurements of the deviation of sin 2/3 e ff 
in AAA modes from sm.2(3j/^, Ks may indicate New 
Physics in b — ► s penguin-induced modes. It is of 



great importance to examine and estimate how much 
of the deviation of sin2/3 fj is allowed in the SM. Ow- 
ing to the presence of color-allowed tree contributions 
in P° — > A + A~As, this mode is subject to a po- 
tentially significant tree pollution and the deviation 
of the mixing- induced CP asymmetry from that mea- 
sured in B — > J/ipKs could be as large as 0(0.10). 
Since the tree amplitude is tied to the nonresonant 
background, it is very important to understand the 
nonresonant contributions in order to have a reliable 
estimate of sin2/3 c ff in AAA modes. 

The deviation of the mixing- induced CP asymme- 
try in P° -> K+K-Ks, K S K S K S , K s n+Tr' and 
KgTr°Tr° from that measured in B — » <ft c£ Kg, i.e. 
sm2(3^ Ks = 0.681 ± 0.025 0, namely, Asin2/3 off = 
sin 2/3 c ff — sin 2f3^ a -K s , is shown in Table IVTIl Our cal- 
culation indicates the deviation of the mixing-induced 
CP asymmetry in B — > K + K~ K$ from that mea- 
sured in B — > cf> cE Ks is very similar to that of the 
A5A5AS mode as the tree pollution effect in the for- 
mer is somewhat washed out. Nevertheless, direct CP 
asymmetry of the former, being of order —4%, is more 
prominent than the latter. 



VII. CONCLUSIONS 

It is important to account for the large nonresonant 
amplitudes in the study of charmless 3-body baryonic 
B decays. We have identified two distinct sources of 
nonresonant contributions: a small contribution from 
the tree transition and a large source of the nonreso- 
nant signal in the matrix elements of scalar densities, 
e.g. (AA|ss|0), induced from the penguin transition. 
This explains the dominance of the nonresonant back- 
ground in B — > AAA decays, the sizable nonresonant 
fraction in A~7r + 7r~ and A 7r + 7r~ modes and the 
smallness of nonresonant rates in B — > tttttt decays. 
The seemingly huge discrepancy between BaBar and 
Belle for the nonresonant contribution in the decay 
B~ — » A~7r + 7r~ is now relieved. Since penguin con- 
tributions to charm decays are GIM suppressed, hence 
nonresonant signals in D decays are always small. 

We have computed the resonant and nonresonant 
contributions to charmless 3-body decays and deter- 
mined the rates for the quasi-two-body decays B —> 
VP and B — > SP. Time-dependent CP asymmetries 
sin2/3 c ff and Acp in A + K~ As, A5A5A5, Kstt + tt~ 
and As7T°7r modes are estimated. Since we have a re- 
alistic model for resonant and nonresonant contribu- 
tions, our estimation of sin 2(3 c g for 3-body B decays 
should be more reliable and trustworthy. 
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TABLE VI: Branching ratios of quasi-two-body decays B — * VP and B — > SP obtained from the studies of three- 
body decays based on the factorization approach Theoretical uncertainties have been added in quadrature. QCD 
factorization predictions taken from [23j] for VP modes and from [24|] for SP channels are shown here for comparison. 
The assumption of £>(/o(980) — > n + n~) — 0.50 has been made for the QCDF calculation. Experimental results are 
taken from [^. Note that the BaBar results for Kq° (1430)7r° and Kg~ (1430)tt + are obtained by neglecting nonresonant 
contributions to (Kty)qtt [l8|] and hence may not be appropriate to compare with Belle directly. 



Decay mode 




BaBar 


Belle 


QCDF 


Theory [111 


<j,K° 




8.4t^g ± 0.5 


9.0±i g ± 0.7 


a 1 +0.4+1.7+1.8+10.6 

— 0.4—1.6 — 1.9— 3.0 


^"^—0.9 


(j)K~ 




8.4 ± 0.7 ± 0.7 


9.60 ± 0.92l O g4 


4 ir+0.5+1.8+1.9+11.8 
a — 0.4— 1.7— 2.1— 3.3 


5.9+ii 

— 1.0 


K"\- 




10.8 ± 0.6li 3 


9.7 ± 0.6 + n a 

— 0.9 


o c+0. 4+1. 5+1. 2+7.7 
—0.3 — 1.4—1.2—2.3 


4.4+" 

— 1.0 


K*°n° 




3.6 ± 0.7 ± 0.4 


< 3.5 


r, 7 +0. 1+0. 5+0.3+2.6 
u - ' — 0.1 — 0.4— 0.3 — 0.5 


1.5+S-S 

^ — . 4 


K*~n + 




II 7 + li 

III -1.2 


8 4±1 


o o+l. 4+1. 3+0. 8+6. 2 
1.2 — 1.2— 0.8 — 1.6 


3 l +0 n 9 o 

°' -0.9 


K*-tt° 




6.9 ± 2.0 ± 1.3 




o o + l. 1+1. 0+0. 6+4.4 
'^" ,3 — 1.0 — 0.9— 0.6 — 1.4 


2.2±° 6 5 


K*°K~ 




< 1.1 




Q on+0. 11+0. 12+0. 09+0. 57 
—0.09 — 0. 10— 0.09 — 0. 19 


0.35+Hr 

— 0.0b 


p°K- 




3.56 ± 0.45 + n 'aI 


3.89 ± 0A7toH 


O ft+0. 9+3. 1+0.8+4.3 
—0.9—1.4—0.6 — 1.2 


— 0.7 


p°K° 




4.9 ± 0.8 ± 0.9 


6.1 ± 1.0 ± 1.1 


A a+0. 5+4.0+0. 7+6.1 
u — 0.5— 2.1 — 0.7— 2.1 


2 o+i? 

^ ,,J — 0.9 


p+R- 




8.0t° s 3 ± 0.6 


ic -1+3.4+2.4 
—3.3—2.6 


7 4+1.8+7.1+1.2+10.7 
' — 1.9— 3.6 — 1.1— 3.5 


— 1.4 


p-K° 




8.0±^3 ± 0.6 




c o+0. 6+7. 0+1. 5+10. 3 
°-° — 0.6— 3.3 — 1.3— 3.2 


1 3+^-S 
i,J — 0.9 


r 




8.8 ± S 

— 0.9 


8.0 + ?,' 3 , ± 0.7 
— 2.0 


i i q+6. 3+3. 6+2. 5+1. 3 
y — 5.0 — 3.1 — 1.2 — 1.1 


7.7+K 

1 ■ 1 —1.6 


r " 








91 9+10.3+8.7+1.3+2.0 
— 8.4—7.2 — 2.3—1.6 


w 3.5 


0~^7V~ 
r " 








1 c 4+8.0+5.5+0.7+1.9 
± °-^ t — 6.4— 4.7— 1.3 — 1.3 


8 5+h 1 , 

1.0 






1.4 ± 0.6 ± 0.3 


3.0 ± 0.5 ± 0.7 


n 4+0.2+0.2+0.9+0.5 
u- — 0.2— 0.1— 0.3— 0.3 


^• KJ — 0.2 


fo(980)K°;fo -< 


7T~'"7r 


5.5 ±0.7 ±0.6 


7.6 ± 1.7 + no 


R 7 +0.1+2.1+2.3 

\J » 1 n "i "i r 1 "i 

— 0.1 — l.O— 1.1 


7.7+S-i 

— 0. * 


/o(980)Jf-;/ - 


* Tr + Tr _ 


9.3±1.0l°j 


8.8 ± Q.8±°{ 9 8 


7 O+0.2+2.3+2.7 
( -O no 1 jR 10 


7.7±g;3 


/o(980)if°;/o -. 


K+K- 


5.3 ±2.2 






c: o+O.l 
O-»-0.5 


/o(980)A-;/ - 


* K+K~ 


6.5 ±2.5 ±1.6 


< 2.9 




7 0+ ' 4 


/o(980)7r-;/ -» 


TT + n~ 


< 3.0 




n c+o. 0+0. 2+0.1 
u - J -o. 0-0. 1-0.0 


U.O»_ Q2 


/o(980)Tr-;/ -» 


K+K~ 








W+ 06 

U.OU_ 04 


/o(980)7r°;/ ^ 


7T + 7T~ 






n n9+°- 01+0. 02+0. 04 

u - uz -o. 01-0. 00-0. 01 


nio+ 003 

U.U1U_ 002 


A' t *°(1430)7r- 




32.0 ± 1.2 +1 °'* 


51.6±1.7±7.4 


-, -, n +10. 3+7.5+49. 9 
11. U_ 6.0-3.5-10.1 


16.91^ 


Fo°(1430)7r° 




13.1±i:|±?;S ±3.6 


9.8 ±2.5 ±0.9 


fi 4+5.4+2.2+26.1 
"■^-3.3-2.1- 5.7 


6.8 + _U 


if*"(1430)7r+ 




28.2±H± 4 3 l ± 5.2 


49.7 ± 3.8±6° 


1 1 q+9.4+3.7+45.8 

^••J-s.s-s.t- 9.9 


i6.2_ 4 ; 


Kl~ (1430)7r° 








+4.7+1.6+22.3 
"J-J_2.8-l.7- 4.7 


8.9l 2 2 ;« 


K *°(1430)i<:- 




< 2.2 






1 0+0.3 



TABLE VII: Mixing-induced and direct CP asymmetries for various charmless 3-body B decays. Experimental results 
are taken from 0|. 



Decay 


sin 2(3 c ft 


A sin 2/3 c ff 


Expt 


A f (%) 


Expt 


K+K'Ks 

KsK s Ks 

K s ty°ty" 


n 798+°' 001+0. 002+0. 009 
U. ( -O_ .002-0 .001-0 .020 
71 Q+0. 000+0. 000+0. 008 
U- 1 l»_o. 000-0. 000-0. 019 
n 79Q+0. 000+0. 001+0. 009 
U. 1 ^»_o. 000-0. 001-0. 020 
fl 71 Q+0.001+0. 017+0. 008 
U. 1 lO-o. 001-0. 007-0. 018 


041 + - 028 

U.U^1_ 033 
u.uoy_ 32 
0.049_S: o£ 
0.038ir o 3 32 


0.05 ±0.11 
-0.10 ±0.20 
-1.20 ± 0.41 


A eq+1.38+0. 53+0.40 
1-UO_ !.oi_o. 54-0. 34 
fl RQ+0-01+0.01+0.05 
U.U»_ .01-0.03-0.07 
n 90+0. 09+0. 07+0. 02 
u - zo -0. 06-0. 06-0. 02 
A Q/l+0. 03+0. 03+0. 32 
*- y *-0. 02-0. 05-0. 40 


-7±8 
14 ± 15 
-18 ±22 
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